Electromagnetic nonlinear X-waves 
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Nonlinear optical media that are normally dispersive, support a new type of localized (non- 
diffractive and nondispersive) wavepackets that are X-shaped in space and time and have slower 
than exponential decay. High-intensity X-waves, unlike linear ones, can be formed spontaneously 
through a trigger mechanism of conical emission, thus playing an important role in experiments. 



The nonlinear response of condensed matter can com- 
pensate for the diffractive spreading of optical beams, or 
the dispersive broadening of pulses due to group- velocity 
dispersion (GVD), forming spatial or temporal soli- 
tons 1^, respectively. Recent experimental results con- 
cerning the self-focusing behavior of intense ultrashort 
pulses 1^, @, H, |[ 0, H] indicate, however, that the spa- 
tial and temporal degrees of freedom cannot be treated 
separately. When the three length scales naturally asso- 
ciated with diffraction, GVD, and nonlinearity become 
comparable, the most intriguing consequence of space- 
time coupling is the possibility to form a nondiffractive 
and nondispersive localized wavepacket (LWP), namely 
a spatio-temporal soliton or light bullet |^ characterized 
by exponentially decaying tails. A strict constraint for 
the formation of light bullets is that the nonlinear phase 
changes counteract both the linear wavefront curvature 
and the GVD-induced chirp, leading to space-time focus- 
ing, as occuring in Kerr-like focusing media with anoma- 
lous GVD g §. 

Viceversa, a normal GVD rules out the possibility to 
achieve bullet-type LWPs. In this regime, the field evo- 
lution is known to be qualitatively different, and involves 
complex phenomena such as temporal splitting and spec- 
tral breaking (H ^, 0] . This is the reason why no attempts 
have been made to answer the fundamental question as 
to whether any form of nonlinearity-induced localization 
could still take place in normally dispersive media. In 
this letter, we show that LWPs do exist also with normal 
GVD in the form of nonlinear X-waves (NLXWs) or X- 
wave solitons. To date X-shaped waves are known only 
in the context of linear acoustic [0 , or electromagnetic 
propagation, and constitute the polichromatic gener- 
alization of diffraction- free Bessel (or Durnin [12|) beams. 
They have been observed in both acoustical 1 13 1 and op- 
tical experiments, both requiring beam-shaping tech- 
niques. Here, we find propagation-invariant NLXWs that 
can be naturally regarded as the continuation of linear 
X-waves into the nonlinear regime. Yet, we find one 
fundamental difference between linear and NLXWs. At 



high intensity, the formation of X-shaped LWPs occurs 
spontaneously from conventional bell-shaped (in space 
and time) beams throughself-induced spectral reshaping 
triggered by a mechanism of conic al emission. Thus, 
NLXWs are expected to have stronger impact on experi- 
ments than linear X waves. 

To support the generality of NLXW concept, we 
choose two different phenomena, whose spatio-temporal 
dynamics have been widely investigated experimentally. 
Specifically, we consider self-action of a scalar wavepacket 
Ml (carrier luq) due to a pure (cubic) focusing Kerr ef- 
fect, or generation of a U2 wavepacket at second-harmonic 
(SH, 2wo) in non-centrosymmetric (quadratic) media. In 
the paraxial regime, the evolution in Kerr media is ruled 
by the scalar 1-1-3 nonlinear Schrodinger (NLS) equation. 



idi^ui + V'^ui — didrrUi + rjuipMi = 0, 



(1) 



whereas SH generation is ruled by the vector NLS model 
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In Eqs. ([T]^) the link with real- world variables AT, F, Z, T 
is as follows: C = -^/^d/ is the propagation distance in 
units of diffraction length Z^f — IkiW^ associated with 
the beam waist VFq, = is the transverse Lapla- 

cian where (AT, F) = Wo(C,??), and t ^ {T - Z/Vgi)/To 
is time in a frame traveling at group- velocity Vgi (of ui) 
in units of To = {\k'(\Zdf 12)^/^ , k'^ (m = 1,2) being 
the GVD at mcjQ. The coefficients are dm = k'^/\k'l\, 
fm — ki/km (ct2 — 1/2), the wavevector mismatch Sk = 
{k2~2ki)Zdf, and F = Z^f /Zni (not rescaled out to quan- 
tify the impact of nonlinearities) . Here Zni — {x3^p)~^ 

are nonlinear 



in Eq. M) and Z„ 



in Eqs. 



length scales associated with the input peak intensity Ip 
at ujo (|ui|^„^ = 1), and X3 [m/W] and X2 [W^i/^j are 
standard nonlinear coefficients. Finally v = ZdfSV/To 
accounts for the walk-off due to group- velocity mismatch 
(GVM) SV = Vg2^ — Vgi^ , though we analyse (where not 
stated otherwise) the GVM-matched case v — 0. 
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We seek for propagation- invariant, radially-symmetric 
LWPs of the form ui = fi{r, t) exp(— i/3C), accompanied 
in Eqs. (||) by a symbiotic SH U2 = f2{r,t) exp[— j(2/3 + 
Sk)(]. Here /i_2 are real, and /3 is a nonlinear phase-shift 
(in Eqs. (^) it makes the two waves nonlinear ly phase- 
matched). In the cubic case, /i obeys the equation 



/i 



7i-^i9«/i + &/i+7/i =0, 



while Eqs. yield (after setting /i/^/ctJ ^ /i), 

/i + r- Vi - Didttfi + bfi + 7/2/1 = 0, 
h + r-'h ~ D2duf2 + a/2 + 7^ = 0' 



(3) 



(4) 



where we have set = c^m/cm, / = df/dr, and 

r = \p\p ^ (x^W?'^ {P^ ^ eW).t= |/3|r, 7 - r/|/3|, 
h — and a — [2b-\-5k/\(}\)/a2 in the nondegenerate 

case {(i ^ 0), while (r, = (p, t), a = 5k/ <72, 7 = F, 
and 6 = in the degenerate case (/3 = 0). Eqs. (||- 
^ must be integrated along with the boundary condi- 
tions /i,2(r-,±oo) = /i,2(oo,<) = 0, and /i,2(0,<) = 0. 
While the anomalous GVD regime {Dm < 0) guaran- 
tees that, for /? < (&, a < 0), nearly separable LWPs 
(light bullets) exist with exponentially decaying tails, in 
the normal GVD regime {Dm > 0) the nature of the 
LWP solutions (if any) must change dramatically, be- 
cause the low-intensity exponential damping no longer 
takes place. Pseudo-spectral numerical techniques (i.e., 
solve Eqs. (^^) as a dynamical evolution problem in r 
with appropriate discretization in t |l5|l ) are well suited to 
search for strongly nonseparable objects with slow spatio- 
temporal decay. We have implemented such methods 
(using two different algorithms), and found LWPs when 
/3 > (&, a > 0). Efficient convergence occurs by em- 
ploying as a trial function (of time, at r = 0) the real 
part of the waveform, 

= ;m = l,2, (5) 



y/{A~t{T/dm.W+e+V 

which represent X-shaped LWP solutions of Eqs. (|l[^) in 
the linear limit (F = 0). Here A is a free parameter: the 
smaller A, the stronger the localization. The LWP solu- 
tions, obtained from Eq. (^) with A = 1 in Eq. (^) are 
shown in Fig. |^. For 6 = and moderate nonlinearities 
[7 = 1, Fig. 1(a)], the "ground-state" LWP mode has a 
clear X-shape (in x — t, or V-shape m r — t variables), 
encompassing a space-time tightly confined structure, 
with slow axisymmetric spatial decay {'^ l/r) accom- 
panied by (radially increasing) temporal pulse-splitting. 
In the case 6=1, while the field maintains its basic X- 
shape, it develops radial (damped) oscillations as shown 
in Fig. 1(b), e.g. for the strong nonlinear case (7 = 10). 
From Eqs. (^ we obtain similar NLXWs, where the two 
symbiotic LWPs /1-/2 can be both of ground-state type 
{b = a = 0), ground-oscillatory type (6 = 0, a ^ 0), or 



FIG. 1: Spatio-temporal field profile 111(2;, y — 0,t) of NLXWs 
in normally dispersive {Di = 1) Kerr media, as obtained from 
Eq. (3) with: (a) 6 = 0, 7 = 1; (b) 6 = 1, 7 = 10. 



FIG. 2: NLXW spatio-temporal intensity (a) \ui{t,r)\'^; (b) 
\u2{t,r)\^, obtained from Eqs. (4) with 6 = 1, a = 2, 7 = 1. 



both of oscillatory type {b = 1, a 7^ 0, with out-of-phase 
coherent oscillation) , an example of the latter case being 
shown in Fig. |^. 

In both nonlinear processes the oscillations stem from 
the fact that the spatial behavior of the low-intensity 
portion of the LWPs is governed by a zero-th order 
Bessel equation, which can be easily obtained by Fourier- 
transforming the linear (7 = 0) limit of Eq. (||) or (de- 
coupled) Eqs. (^. Therefore NLXWs exhibit the char- 
acteristic damped spatial oscillations of Jq Bessel func- 
tions, and can be regarded as the natural generalization 
of monochromatic nondiffractive conical Jq beams, whose 
physical realization was demonstrated both in the linear 
|]l2| and nonlinear [|6| regimes. More precisely, once the 
NLXW solutions /„, = fm{r, t) are obtained, they can be 
represented in the spectral domain of transverse wavevec- 
tor K = {K^+Ky Y^^ and frequency detuning fi, through 
the Fourier-Bessel transform Sm = Sm(K, il), as 



/m — 



-f 00 



+00 



Sm{^,K)Ja{Kr)e'^^*KdKdn, (6) 



showing their nature of a weighted superposition of co- 
nical Jo beams with different frequency. 

Compared with their counterpart in the anomalous 
GVD regime (i.e., 1-H3D bullet LWP), NLXWs possess 
important differences: (i) NLXWs experience weaker lo- 
calization, i.e. l/r instead of exponential decay. As 
a consequence the energy of X-waves (||) and likely 
also of our NLXW solutions, is infinite, viz Ei = 
Jq rdr J^^dt\ui\^ in Eq. (||) tends to infinity as the 

computational window a ^ 00; (ii) unlike bullets that 
vanish in the linear limit, NLXWs have a finite limit for 
F = 0; (iii) NLXWs are not unique, in the sense that 
inifinitely many solutions can be found for fixed P and 
F. This follows also from (ii) since also linear paraxial X- 
waves are a continuous family, e.g. parametrized by A in 
Eq. (H); (iv) NLXWs and bullets exist for different sign 
of (3, entailing opposite phase shifts in Kerr media, and 
different constraints {(3 < -~Sk/2 for NLXWs, (3 < 5k/2 
for bullets) in Eqs. (^. 

Once the existence of nonlinear LWPs supported by 
normal GVD is established, one might wonder about 
their importance and observability. In the linear regime, 
both Durnin Jq beams [|l^ and nonmonochromatic X- 
shaped LWPs |l^ can be observed only by means of 
experimental arrangements that adapt the input to the 
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FIG. 3: Comparison of NLXW and CE gain in the spectral 
domain: (a) \Si{K,^l)\ transform of the field uj in Fig. 2(a) 
{\S2(K,Q,)\ is similar); (b) CE domains for the out-of-phase 
cw plane- wave eigenmode with the same parameters. 



FIG. 4: Output intensity profile |-Eip = /p|itip, as obtained 
from Eqs. (2), (see text for coefficient values) in the GVM- 
matched case n = (top panel), and with GVM, v ^ (bot- 
tom panel). 



LWP (e.g., by means of lensacons or axicons). Vice- 
versa, we have found that the interplay of the non- 
linearity and normal GVD is responsible for a univer- 
sal mechanism, namely colored conical emission (CE), 
that allows for the self-induced spectral (in fl-K) re- 
shaping necessary to turn conventional, e.g. gaussian, 
pulsed beams into X-waves. In fact, it is known that 
in Kerr media [0 the (modulational) stability analy- 
sis of the cw plane-wave solution u = e*'"^ of Eq. (1) 
yields exponential amplification of conical plane-wave 
(or Bessel) perturbations with wavevector K and fre- 
quency detuning fl, such to yield real values of gain 
g = {[{Din"^ - K'^ f + Tf - r2}i/2. This analysis can be 
readily generalized for the cw plane-wave eigensolutions 
Ml = uio exp(— «/3C) U2 = U20 exp[— z(/3-|-(5fc)C] of Eqs. (2- 
3) H to obtain the gain g^[h±{b'^ - c)l/2]l/^ where 
2h = T{uIq - 2ulo) - ^1 - ^2, c ^ {^1^2 - ru2g)2 - 
Vu\q^\, and fim = ^mK'^ — dm^l"^ + m(3 — (m — l)6k, 
m = 1,2. In the normal GVD regime {dm > 0) these 
expressions entail CE, i.e. preferential amplification of 
waves at K (angles) linearly increasing with frequency 
detuning. By comparing (see Fig. ||) the CE gain with 
the NLXW spectrum, obtained by inverting Eq. (^) , it is 
clear that the instability provides amplification at K ^fl 
pairs that favour the formation of NLXWs. Although the 
stability analysis is carried out for cw plane-wave pum- 
ping beams, CE occurs also from tightly-focused short- 
pulse input beams |l7|. In this case, we expect CE to 
amplify frequencies K — il contained in the broad in- 
put spectrum, while preserving the phase coherence be- 
tween different spectral components necessary for the 
formation of a NLXW. In order to support this con- 
jecture and prove that NLXWs are the key to under- 
stand the dynamics of experiments carried out with nar- 
row beams and short pulses in the normal GVD regime, 
we have conducted numerical simulations of the propa- 
gation. This is crucial also to assess the observability 
of such type of LWPs with finite-energy, which are the 
nonlinear counterpart of monochromatic finite-aperture 
Jo or Bessel-Gauss beams, and nonmonochromatic 
finite-energy linear X-waves observable in the 

real world. While extensive results will be reported else- 
where, we focus specifically on quadratic media [Eqs. (|^)] 
in the large negative mismatch limit, where the field ui 
plays a leading role, and experiences an effective focus- 
ing Kerr effect |^. This case has twofold relevance: 
(i) higher-order effects not included in our models (Ra- 
man, self-steepening, space-time coupling, saturation, 
etc.) have lesser impact on the propagation as com- 
pared with true Kerr media; (ii) experimental results in 



Lithium Triborate (LBO) indicates the occurence of pulse 
compression in spite of the fact that the medium is nor- 
mally dispersive p9| . We model the latter case assuming 
X2^7x 10-5 k'{ ^ 0.015 psVm, k'^ = 0.07 psVm, 

a mismatch Afc — —30 cm"^ {5k = —180), and a spatio- 
temporal input gaussian beam ui{Z = 0) with FWHM 
170 fsec duration and 65 fim beam width. Figure || (top 
panel) shows the output intensity in a 4 cm crystal, in 
the GVM-matched case (A = 1.3//m, T = 35 C in LBO) 
for an intensity Ip — 50 GW/cm^ (Z„; — 0.6 mm). As 
shown, while a moderate fraction of the energy lags be- 
hind and in front in the form of pulse satellites, the main 
portion of the beam develops the characteristic structure 
of a NLXW. Importantly, we observe that the forma- 
tion of NLXW-type of LWPs in Fig. ^ is accompanied by 
strong pulse compression. In fact, it can be shown that 
the peculiar spatio-temporal structure of a NLXW leads 
to an effective GVD, which we obtain in the form A:"^^ = 

^ ^ SlH^) - 4] ^ ^" - fco(£)'' ^liere the angle 
9 ~ sin 6* = kt/ko, and kt = (fcg - k^Y^^ = K/{\l3\Wo). 
Therefore the dispersive contribution which stems from 
the angular dispersion 9 — 9{uj) counteracts the material 
(normal) GVD, leading to an effective anomalous GVD, 
which in turn explains the compression. 

NLXWs are robust also against strong GVM. First, 
when w 7^ in Eqs. (|), we find CE as well as NLXW 
solutions of Eqs. (||) characterized by an additional com- 
plex phase profile (details will be given elsewhere). On 
the other hand, our simulations show that, at sufhciently 
high intensity, after a short distance, the launched (ui) 
and generated (U2) fields tend to develop spontaneously 
NLXW shapes, meanwhile leading to nonlinear walk-off 
compensation (ui,2 travel locked together). For instance 
in Fig. ^ (bottom panel) we show the profile ob- 
tained for an input intensity Ip = 70 GW/cm^ after 1.5 
mm propagation in LBO at A = 1.06/J,m, where a GVM 
as large as 5V = 45 ps/m {v = 75) is compensated. 

Our results show that, contrary to common belief, non- 
linear space-time localization takes place also in normally 
dispersive media. NLXWs are the eigenmodes of 1-I-3D 
paraxial wave propagation models, and will be central to 
interpret correctly numerical and experimental results. 
Having restricted ourselves to NLXWs that travel at the 
natural group- velocity of light, further work will be de- 
voted to search for NLXW solutions with finite energy (of 
which we have given numerical evidence), and/or sub- or 
super-luminal nature, as well as the role played by non- 
paraxiality and higher-order nonlinear terms. 
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